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Abstract 

New Faddeev-Yakubovsky calculations are reported for AA He and /^Be in 
terms of clusters of a's and A's, using AA s-wave potentials motivated by 
several of the Nijmegen model interactions. The self consistency of the AA 
hypernuclear world data for these species is discussed. The newly reported 
Aj ^He event is found to be compatible with AA interaction strengths provided 
by the Nijmegen soft-core one-boson-exchange model NSC97. Faddeev cal- 
culations for AA H and AA He suggest that these AA hypernuclei are stable 
against emitting A's for any (essentially attractive) AA interaction, whereas 
calcualtions for AA H do not allow a clear-cut conclusion whether or not it is 
particle stable. 
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I. INTRODUCTION 



Data on strangeness S = — 2 hypernuclear systems is rather scarce, and no data exist 
for systems with higher strangeness content of hyperons (Y). Multistrange hadronic matter 
in finite systems and in bulk is predicted on general grounds to be stable, up to strangeness 
violating weak decays (see Ref. |l|] for a recent review). Until recently only three candidates 
existed for AA hypernuclei which fit events seen in emulsion experiments |2|-Q. The AA 
binding energies deduced from these events indicated that the AA interaction is strongly 
attractive in the 1 S channel |5]||, although it had been realized J7]|| that the binding 
energies of the two older events, ^He and A^-^e, are inconsistent with each other. This 
outlook is perhaps undergoing an important change following the very recent report from 
the KEK hybrid-emulsion experiment E373 on a new candidate U for AA He, with binding 
energy substantially lower than that deduced from the older event ||. Furthermore, there are 
also indications from the AGS experiment E906 for the production of light AA hypernuclei 
|T0| , perhaps as light even as AA H, in the (K~ , K + ) reaction on 9 Be. 

The study of multistrange systems can provide rather stringent tests of microscopic mod- 
els for the baryon-baryon interaction. Over the years, the Nijmegen group has constructed 
a number of one-boson-exchange (OBE) models for the baryon-baryon interaction (NN, 
AN — UN, and EN — AA — EE — AS) using SU(3)-flavor symmetry to relate baryon-baryon- 
meson coupling constants and phenomenological hard or soft cores at short distances (for a 
recent review, see Ref. |]TT|). In addition, the Jiilich group has constructed OBE models for 
the YN interaction along the lines of the Bonn Model for the NN interaction, using SU(6) 
symmetry to relate coupling constants and short-range form factors |L2|. The SU(6) quark 
model has also been used to derive baryon-baryon interactions within a unified framework 
of a (3g) — (3g) resonating group method, augmented by a few effective meson exchange 
potentials of scalar and pseudoscalar meson nonets directly coupled to quarks |H| . In all of 
these rather different baryon-baryon interaction models only 35 YN low-energy, generally 
imprecise data points serve for the purpose of steering phenomenologically the extrapolation 
from the NN sector, which relies on thousands of data points, into the strange YN and 
YY sectors. It is therefore of utmost importance to confront these models with the new 
AA hypernuclear data in order to provide meaningful constraints on the extrapolation to 
S = — 2 and beyond. 

In this paper we report on Faddeev-Yakubovsky three-body and four-body calculations of 
light AA hypernuclei, using AA s-wave potentials which simulate the low-energy s-wave scat- 
tering parameters produced by several Nijmegen OBE models. Our own work is compared 
to other recent variational calculations [p^ -[T6||. The advantage of the Faddeev-Yakubovsky 
method is that it accounts properly for all the relevant rearrangement channels. The pur- 
pose of our calculations is twofold: to check the self consistency of the data, particularly 
for AA He and AA Be |T7| which are treated here as clusters of a's and A's; and to find out 



which of the recent Nijmegen Soft-Core (NSC97) models JT8|,|19j, or the Extended Soft-Core 



(ESCOO) interaction model [ [Ll[| , is the most appropriate one for describing well AA hypernu- 
clei. The paper is organized as follows. The Faddeev-Yakubovsky approach, its restriction 
to s waves for light hypernuclei, and some computational aspects are described briefly in 
Sect. [TT| We specify the input potentials, fitting as much as possible known binding energies 
of given hypernuclear subsystems. The results of our Faddeev-Yakubovsky calculations for 
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Ai \He and \^Be in terms of clusters of a's and A's, using AA s-wave potentials motivated by 
several Nijmegen model interactions, are reported and discussed in Sect. |TTl[ Faddeev calcu- 
lations for AA *H and A^He are also reported. A preliminary, partial report of these results has 
just appeared as a Rapid Communication ||20|| . In the Summary section we briefly discuss 
the implications of our work on the stability of aa-H, for which our Faddeev-Yakubovsky 
calculations at present are in progress. 



II. METHODOLOGY 

A. Faddeev equations 

The bound states of the three-body systems considered in this work are obtained by 
solving the differential Faddeev equations 
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(H + Vp c + V^-E)Up = -V^ £ U, , 0=1,2,3 , (1) 

where and V* are Coulomb and short-range pairwise interactions respectively in the 
channel (3, Ho is the kinetic energy operator, E is the total energy and the wavefunction of 
the three-body system is given as a sum over the three Faddeev components, \I/ = J2p=i Up, 
corresponding to the two-body breakup channals. The Faddeev components are functions 
of spin-isospin variables and of relative Jacobi coordinates. When two constituents of the 
three-body system are identical, for example A hyperons in A^He (aAA) generically of the 
form CAA (C = core), the coupled set of Faddeev equations simplifies as follows: 

(H + V AA - E)U = -V AA (W - P 23 W) } (H + V AC - E)W = -V AC (U - P 23 W), (2) 

where P23 is the permutation operator for the A hyperons (particles 2,3). U is the Faddeev 
component corresponding to the rearrangement channel a — (AA) and W corresponds to the 
rearrangement channel (ah) — A: 

*a-(AA) = U , * (aA )-A = (1 - P2Z)W ; $ = ^ a _( A A) + *(oA)-A • (3) 

The probability weights for the various rearrangement channels (3 in the total wavefunction 
^ of ^ a He are defined by 

< ^ > 

^TW' (4) 

For ^Be, here considered as a Acta system, a similar simplification occurs, but where the 
minus sign in Eqs. (§) in front of P23 arising from the Fermi-Dirac statistics for the A's is 
replaced by a plus sign owing to the Bose-Einstein statistics for the a's, and V c is added 
on the l.h.s. of the first Faddeev equation. These Faddeev equations, for ^Be(| + ) and 



AA J 



He(0 + ) ground states, were solved in Ref. using an s-wave approximation in which 



only the £ a = 0, \ a = partial waves were retained for the Jacobi-coordinate x a , y a degrees 
of freedom, respectively, where 
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x a = r /3 -r 7 , y a = ■ r a , (3 > 7 ^ a, a, (3, 7 = 1, 2, 3, (5) 

m^ + m 7 

where m a and r a , a=l,2,3, are the mass and the position vector of the a-th constituent, 
respectively. (The total orbital angular momentum is given here by L = £ a + X a ). The 
present calculation extends and updates the previous one. We comment in Sections |TTTj and 
[TV] on the validity of the s-wave approximation. 

Other three-body systems studied in the present work are the isodoublet aa^, AA^e 
charge-symmetric hypernuclei, here considered as 3 HAA and 3 HeAA respectively. For the 
ground state (| + ) of this system, the s-wave Faddeev equations assume the form 

1 

f *^i2/i 

(Ho{xi,yi) +Vaa(zi) -E)U(xi,yi) = -VaaC^i) / ^itt^^UI). 



{K(x 2 ,y 2 ) + V AC (x 2 ) - E)W(x 2 ,y 2 ) = -^V AC (x 2 ) 



1 

dv 2 ^A T U(x' 2 ,y' 2 ) (6) 

x 2 y 2 



1 

+ [ dv 2 ^BW(x>>,y») 
J x%m 



x 2 y 2 

where v s Kk (x) is the singlet AA potential, Vac = diag{v AC , v AC } with v s AC and v AC the singlet 
and triplet AC interaction potentials, respectively, 

*=c4-f>- 

W(x,y) = {W{x,y),W(x,y)) T . 
The appropriate reduced masses are given by 

m c m A m A (m c + m A ) m A m c m A 

M l = 2 n; > M 2 = — " , /Xl = — , A*2 " 



mc + 2mA "^c + 2mA 2 mc + mA 

The appropriate transformation of coordinates is given by 

/ // A^ \2 1 2 , o ^ \l/2 / m C //^A N 2 . 2 n mA \l/2 

^1 = — xx) +yi +2 — xxyxvx) 1 , y x = ■ ((tt x V + ^ ~ V^ 1 ^ 1 ^ , 

mA mA m c + m\ Mi Mi 

11 // A 4 ! \ 2 1 2 o A* 1 \l/2 // m C f/f^K \ 2 1 2 1 mA \l/2 

»i = xi) +Vi - 2 x 1 y 1 v 1 ) 1 , yj = ■ — zi +Vi + 2—x 1 y 1 v 1 ) 1 , 

mA mA mc + mA Alj M± 

1 11 A t 2 \ 2 1 2 1 o A*2 \l/2 / l// m A A 2 . 2 q m A U/2 

^2 = ^2 + 2/2 +2 X2I/2V2 7 , V 2 = o T7- X2 - 2 XF a; 2y2W2 ' , 

mA mA 2 M2 M 2 

11 (I A*2 x2 , 2 o A t 2 sl/2 // m A // m C n 2 . 2 , r, 77 ^ M/2 

^2 = £2 +2/2 -2 x 2 y 2 v 2 ) 1 , |/ 2 = ■ — x 2 +2/2 + 2— x 2 2/2^2 ' ■ 

m c m c m c + m A M 2 M 2 
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Here Vj is the cosine of angle between xj and yj. The s-wave total wavefunction then assumes 
the form 

^(x 1 ,y 1 ,v 1 ) =U(x 1 ,y 1 )/(x 1 y 1 )+AW(x[,y[)/(x[y[)+AW(x'l,y'l)/(x';y'l) , (7) 

where the product AW ensures that the two A's are in a 1 S'o state as required by the Pauli 
principle. Note that the squares of elements of A correspond then to a (2J+ 1) average over 
the J = 0, 1 (singlet and triplet, respectively) \H and ^He states. We also note in passing 
that the s-wave Faddeev equations for AA B.e (aAA) and for 9 A Be (Aaa) are obtained from 
Eqs. (P) by setting A = B = 1 with obvious changes in the notation for the masses and the 
total wavefunction (and adding the Coulomb potential between the two a's in A Be). 

B. Yakubovsky equations 

The four-particle wavefunction \1/ is decomposed, in the Faddeev- Yakubovsky method 



21| , into components which are in a one to one correspondence with all chains of partitions. 
The chains consist of two-cluster partitions a 2 {e.g., (ijk)l or (ij)(kl)) and three-cluster 
partitions a 3 (e.g., (ij)kl) obeying the relation 03 C a 2 . The latter means that the partition 
03 can be obtained from the partition a 2 by splitting up one subsystem. Any three-cluster 
partition 03 may be obtained in this way from three different two-cluster partitions 62, C2 
and d 2 . It is easy to see that there exist generally 18 chains of partitions for a four-particle 
system. The components ^ az a 2 obey the differential Yakubovsky equations 



(H + V a , - E)m a3a2 + V a3 ]T # C3(l2 = -V a3 J2 E • ( 8 ) 

Here H is the kinetic-energy operator and V as stands for the two-body potential acting 
within the two-particle subsystem of a partition 03 (e.g., V a3 = Vy with a 3 = (ij)kl). 

We demonstrate the Yakubovsky equations for the case of two pairs, each consisting of 
identical constituents. To be more specific we consider the pair (1,2) to consist of two identi- 
cal fermions, and the pair (3,4) to consist of two identical bosons, altogether simulating the 
four-body AAaa system. This results in a reduction of the number of Faddeev- Yakubovsky 
components into seven independent components satisfying the following Yakubovsky differ- 
ential equations: 

{H + V C + V aa - E)Ui + V aa (I + P U )U 2 = -V aa ((I + P M )U 4 - P 12 (I + P U )W 3 ) , (9) 
(H + V C + V Aa - E)U 2 + Va^U, + P U U 2 ) = -V Aa (P 34 U 4 - P 12 U X + W 2 - P 12 W 3 ) , (10) 
(H + V C + Vaa - E)U 3 + V AA (I - P 12 )U 4 = -V AA ((I - P l2 )U 2 + (I - P 12 )P U W 3 ) , (11) 
(H + V C + V Aa - E)U 4 + V Aa (U 3 - P 12 U,) = -V Aa (P34C/3 - iW 2 + W 1 + P U W 3 ) , (12) 
(H + V C + V AA - E)W 1 + V AA W 2 = -V AA (I - P 12 )U 1: (13) 
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(H Q + V C + V aa - E)W 2 + V aa Wt = -V aa (I + P U )U 3 , (14) 

(H + V C + V Aa - E)W 3 - V Ka P X2 P u W z = -V Aa (U 4 + U 2 ), (15) 

where Up (/3=1,2,3,4) are Yakubovsky components corresponding to three-cluster partitions 
evolving from 3+1 type breakup into two-cluster partitions, W 7 (7=1,2,3) are Yakubovsky 
components corresponding to three-cluster partitions evolving from 2+2 type breakup into 
two-cluster partitions: 

U\ = #(34)2,1 ) U 2 = #(24)3,1 5 Us = #(12)4,3 , ^4 = #(24)1,3 5 (16) 
Wl = #(l 2 ),34 , W 2 = # (3 4),12 , W 3 = # (13 ),24 , (17) 



in the shorthand notation of Ref. [^4] for the (03,02) chains of partitions. For each 
Yakubovsky component the natural Jacobi relative coordinates are employed. For the par- 
ticular AAaa problem under consideration, V c stands for the Coulomb interaction between 
particles 3 and 4, expressed in terms of the coresponding Jacobi coordinates. There are four 
rearrangement channels altogether, defined in terms of the following wavefunctions: 

# (aa A)A = (/ - Pu)Ux + (/ + P34 - P12P34 - Pn)U 2 , (18) 

#(aAA)a = (/ + ^34^3 + (/ + ^34 " ^12^34 " Pu)U 4 , (19) 
#(aa)(AA) =W X + W 2 , #( a A)(aA) = (/ + ^34 " ^12^34 " ^12)^3 • (20) 

The total wavefunction can be written as 

# = #(a«A)A + #(c*AA)a + #(aa)(AA) + #(aA)(aA) • (21) 



The probability weights for the various rearrangement channels /3 in the total wavefunction 

10i 
AA J 



of A^Be are given by Eq. (H) of the previous subsection 



C. Computation 

The s-wave Faddeev- Yakubovsky equations are solved using the cluster reduction method 
(CRM) which was developed by Yakovlev and Filikhin [2lJ and has been recently applied 
to calculate bound states and low-energy scattering for systems of three and four nucleons 
26| . In this method, the Faddeev- Yakubovsky components are decomposed in terms of the 



eigenfunctions of the Hamiltonians of the two- or three-particle subsystems. Due to the 
projection onto elements of an orthogonal basis, one obtains a set of equations with effective 
interactions corresponding to the relative motion of the various clusters. Thus, the Faddeev 
components U a , a=l,2,3 are written in the following form: 

00 

U a (x,y)=Y,<P«,i(x)f*,i(y) ■ (22) 
1=0 
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Here, the basic functions 4> a j are the solutions of the two-body Schrodinger equation for the 
subsystem a: 

h 2 

(-7; 9l + V a (x))(j) a! i(x) = e a j(f) a! i(x) , (f> a ,l{0) = <fia,l( R cutos) = , (23) 

where \x a is the reduced mass of the subsystem a, and v a (x) is the s-wave projected potential 
V a . Substituting (|22]) into ([[]) and projecting onto the orthogonal basis (f) at i(x) one obtains 
a set of integro-differential equations for the functions f a ,i(y)- A fairly small number of 
terms is sufficient, using limited computer resources [PEfl , to generate a stable and precise 
numerical solution independently of the value adopted for -R cu toff, provided the latter was 
chosen sufficiently large. The convergence of a typical Ai\He Faddeev calculation, and of a 
typical AA-Be Yakubovsky calculation, as function of the number of terms N retained in Eq. 
(P^) above is demonstrated in Fig. IB. Such calculations usually require less than 10 terms 
to reach good convergence. The CRM has been recently used by one of us to study \R as 
a three-body pnA system pTfl , and A Be and ^B.e in terms of three-cluster systems Aaa 



and aAA, respectively p2| . Here we update and extend that calculation, while applying 
the CRM for the first time to the solution of the actAA four-body problem in terms of 
Yakubovsky equations. 



D. Potentials 

The AA interaction potentials in the 1 So channel which are used as input to the Faddeev 
equations are of the three-range Gaussian form 

^AA = E^ ) exp(-r 2 /A 2 ) , (24) 



following the work of Hiyama et al. [|1J] where a phase-equivalent AA potential of this soft- 



core form was fitted to the Nijmegen model D (ND) hard-core interaction |28| assuming 



the same hard core for the AW and AA potentials in the 1 5'o channel. Our simulation of 
the ND interaction potential is specified, for 7 = 1, in Table |. For other interactions we 
have renormalized the strength of the medium-range attractive component (i = 2) of this 
potential such that, using values of 7 7^ 1, it yields values for the scattering length and 
for the effective range as close to values prescribed by us. In Table [TI| we list the low- 
energy parameters for several Nijmegen potentials, particularly the soft-core NSC97 model 
AA potentials [HJ, and in parentheses the corresponding values for potentials of the form 
(H]) using the listed values of 7. Three of these fitted 1 Sq potentials are shown in Fig. 0. 
The strongest AA potential is the one with the largest absolute value of scattering length, 
due to the ESC00 model |TJ], and the weakest one is associated with the smallest absolute 



value of scattering length, due to one of the NSC97 model potentials. The difference between 
the various NSC97 potentials {a)-(f) is of second order, all of them being only fairly weakly 
attractive. 

The aa short-range interaction and the Aa interaction are given in terms of two-range 
Gaussian (Isle) potentials 
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V® = V® exp(-r 2 //3 r ( g 2 ) - Kt? exp(-r 2 //^ 2 ) , (25) 



where the superscript j specifies the interacting system. For the aa short-range interaction 
potential we used the s-wave component of the Ali-Bodmer potential |[29|| . A Comment 



about the validity of using this phenomenological potential, which does not reflect at short 
aa distance the underlying inner structure of the a clusters, is made below in Sect. fH\ . 
The Aa interaction potential is due to Ref. |3(J where the binding energy and mesonic weak 
decay of A He were studied. The resulting A-hyperon density distribution has been shown 
very recently |3l| to closely resemble that due to a microscopic calculation of A He using YN 



interactions which simulate those of model NSC97 |18[ . The parameters of these potentials 
are listed in Table |TJ. Similar potentials were constructed by us for the A— 3 H and A— 3 He 
singlet and triplet interactions by refitting the range parameter (3 rcp of the Aa Isle potential 
(p5p to reproduce the known binding energies of A H and \H.e, respectively. These fitted 
range parameters and the binding energies calculated using this procedure are shown in 
Table [TV], in comparison with the experimentally known values • The A Be calculation 

is discussed below in the next section. 



III. RESULTS AND DISCUSSION 



We first applied the aa and Aa potentials, specified by Eq. (|25| ) and in Table [TTlJ to 



the solution of the s-wave Faddeev equations for the aaA system. The calculated ground- 



state binding energy of A Be in this three-body model is given in Table [iy| and is in close 
agreement with the measured value, with no need for renormalization or for introducing 
three-body interactions. Had we used a purely attractive Aa interaction potential such 
as the single- Gaussian form due to Bando et al. [j34j] , the calculated ground-state binding 
energy of A Be would have come out too low by about 1 MeV |2~2|. On the other hand 



the calculations by Motoba et al. |p5| , using an Isle potential (similar to ours) for the Aa 
interaction, overbind A Be by about 0.9 MeV, part of which must arise from going beyond the 
s-wave dominant (£, A) = (0,0) component (98.0% probability weight). This overbinding of 
A Be was subsequently cured phenomenologically by Hiyama et al. 0] by weakening the p- 
wave component of the underlying AA interaction; alternatively it can be cured introducing 
a repulsive Aaa interaction potential to simulate a likely contribution due to ANN three- 
body forces Since our interest in A Be in this work is mostly for using it as a subsystem 
input to calculating AA^e, we should strive to use a model that reproduces as closely as 
possible the binding energies of such subsystems. The significance of this requirement, in 
the framework of a Faddeev calculation for A H, was emphasized in Ref. |36[| . Indeed, our 



s-wave Faddeev calculation for A Be nearly reproduces the observed binding energy. The 
effect of the higher partial waves (£, A) = (2, 2), (4, 4) ignored by us must be largely cancelled 
out in the evaluation of ABA^(^Be) as defined later on. 

Having thus gained confidence in the appropriateness of the aa and Aa input potentials, 
we then applied these potentials to the solution of the Faddeev- Yakubovsky equations for 
A ^He and A^-^e, using different AA interactions generically of the form 
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A. AA 6 He and A ^Be 



The ground-state energies Eaa obtained by solving the s-wave three-body (aAA) Faddeev 
equations for AA He and the s-wave four-body (aaAA) Yakubovsky equations for AA Be are 
given in Table [V| for AA potentials ordered according to their degree of attraction. The 



E\\ for A ^He and B 



AA 



corresponding A A binding energies Baa are given by Baa = 
— (Baa + 0.1 MeV) for ^Be. The strongest AA attraction is provided by the potential of 



the uppermost row, simulating the very recent ESC00 model O] which was partly motivated 
by wishing to get a B A a value for AA^e as close to that reported by Prowse 0; indeed our 
calculation reproduces it. Significantly lower Baa values are obtained for our simulation of 



model ND which, however, reproduces well the Baa value reported for AA Be fll7 |. These 
calculated -Baa values were obtained using a AA potential of the form (p4j), very similar to but 
not identical with the one (also marked ND) used in the cluster calculation |14| which is listed 
separately in the table. The appropriate low-energy parameters, also shown in Table [V], are 
very close to each other, and the variation thereby induced in the calculated B\\( A ^B.e) is 
merely 0.03 MeV within our own s-wave calculation. Furthermore, our calculation essentially 
agrees with that of Ref. [14[] for AA^e, whereas for AA^e our calculation provides about 0.5 
MeV more binding. In view of the restriction of our calculation to s waves, its rough 

13(1 for A^Be should be considered satisfactory. We note that 



agreement with that of Ref. 



14 



for BAAC^Be), 



as function of the 



the variation within the various calculations of Ref. 
microscopically-motivated Aa interactions used by these authors, exceeds 0.5 MeV. We feel 
that typical errors due to limiting our calculations to s waves are a fraction of 0.5 MeV 
for Ai\He (as discussed in Sect. pV|) and about 0.5 MeV for aaBc In addition to the 
simulation of the bare AA ND interaction discussed above, we also listed in Table binding 
energies calculated for a simulation of the G-matrix constructed in Ref. out of the ND 
bare interaction. This was done for the sole purpose of comparing below our calculation to 
those of Refs. [1^,37]; of course, using a G-matrix interaction rather than a bare interaction 
within a full-space calculation is fraught with double counting and should be avoided as 
properly acknowledged in the subsequent work of Ref. [14]. Down the list in Table 0, the 
NSC97 models are seen to give yet lower Baa values, which for AA^e are close to the very 
recent experimental report ||. In between the rows pertaining to model ND and those for 
model NSC97, just for interpolation, we also listed results for two intermediate-strength AA 
interaction potentials specified by their low-energy parameters. 



In Tables [VI] and |VII we list information on the geometrical sizes of AA-He and aa'Bb, 



respectively. As expected, the stronger the AA interaction is - the smaller are the geometrical 
sizes of these few-body systems. We note, without offering any explanation for it, that almost 
all of our calculated sizes, particularly for A^-^e, are significantly larger than those calculated 
in Ref. |14| which is the only comprehensive calculation accounting for all rearrangement 
channels to compare with. An interesting exception is provided by the r.m.s. distance 
between the two a clusters in A^Be: R aa = 3.4 fm in both calculations for model ND, in 
spite of the radically different aa potentials used in these calculations. With such a large 
value of R aa , the two a clusters apparently are little affected by the inner part of their 
potential, which for a more compressed configuration would require a better theoretical 
guidance than provided by the Ali-Bodmer potential |29|] used by us. This insensitivity to 
the inner structure of the a clusters, when A's are added, no longer holds when nucleons are 



S 



added; for example, for the 3a system, where R aa = 3.0 fm in typical C + ground-state 



calculations. Thus, Fedorov and Jensen |38], using the Ali-Bodmer aa potential, had to 



introduce a strongly attractive 3a potential in order to nearly reproduce the energies of the 
two lowest + states in 12 C, whereas Hiyama et al. |L4]], excluding Pauli-forbidden aa states 



from their aa potential, had to introduce a strongly repulsive 3a potential. 



For Ai\He we also listed in Table [V| the probability (f|) of having it as a two-body (a A) — A 
cluster. This probability goes to one as V AA approaches zero; the structure of AA He is then 
well approximated by a A He — A cluster, and the r.m.s. distance R Aa between one of the A 
hyperons and the a cluster is only slightly smaller than the value R Aa ( A He) = 3.16 fm. On 
the other hand, the value of the analogous r.m.s. distance between one of the A hyperons and 
the cm. of the aa clusters in ^Be, R( aa )A = 3.26 fm as shown in Table [VI I| for Vaa —* 0, 
is considerably larger than the value R^a^ijjie) = 2.88 fm, indicating that the two-body 
clustering A Be — A is far from saturating the structure of AA Be in this limit. This point, as 
well as the information given in parentheses in Table |V11| , will be discussed in more detail 
towards the end of the subsection. A final point in discussing Tables [VI| and [VI 1| concerns the 



ground-state energies E AA calculated for the G-matrix interaction ND(G) based on the bare 
ND interaction. Our calculated binding energy for AA Re is higher by about 0.8 MeV than 
that calculated in Refs. [|T(],[J7]]. Given the fact that our calculation for the ND interaction 



comes very close within that of the complete calculation of Hiyama et al. [14], we conclude 



that the calculations of Refs. ||16| , |37| are far from being complete. 

It has been shown in previous cluster calculations 0|| that the calculated Baa values 
for AA He and for AA Be are correlated nearly linearly with each other. Our calculations also 
indicate such a correlation, as demonstrated in Fig. |3| by the dotted line. This line precludes 
any joint theoretical explanation of the AA He an d AA-^e experimental candidates listed in 
Table 0. If ^(^Be) = 17.7 ± 0.4 MeV @[TJ, then the theoretically implied B AA value 
for Ai\He is about 9.1 ± 0.4 MeV, considerably below the value reported by Prowse [f| but 
considerably above the value reported very recently by Takahashi et al. M. If -Baa(a 1 a^ 6 ) = 
14.6 ± 0.4 MeV, on the assumption that the 7r~ weak decay of AA Be ground state occurred 



to the first excited doublet levels of A Be at 3.1 MeV |39], then the theoretically implied 



Baa value for AA He is 6.1 ± 0.4 MeV, no longer in such a spectacular disagreement with 
experiment but still significantly below the value reported by the recent observation |J of 
AA 6 He. 

For Vaa = 0, the lower-left point on the dotted line marked '4-body model' in Fig. |3| 
corresponds to approximately zero incremental binding energy AB AA for AA He, where 

AB AA ( A A A Z) = B AA ( A A A Z) - 2B A ( (A A 1] Z) . (26) 
This is also explicitly listed in Table |V1U| and is anticipated owing to the rigidity of the 



a core. However, the corresponding AB AA value for AA Be is fairly substantial, about 1.5 
MeV, reflecting a basic difference between the four-body aaAA calculation and any three- 
body approximation in terms of a nuclear core and two A's as in AA He. To demonstrate 
this point we show by the dot-dash line marked '3-body model' in Fig. |3] the results of 
a three-body calculation for AA Be in which the 8 Be core is not assigned an aa structure. 
In this calculation, the A— 8 Be potential was chosen in a Woods-Saxon form, with depth 
Vq = 18.12 MeV, radius R = 3.24 fm and diffusivity a = 0.60 fm. The geometry and depth 
were fixed by requiring (i) that the ground-state binding energy of the A— 8 Be system agrees 
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with that calculated in the aaA model calculation for B\(jJ5e) (see Table |V|); and (ii) 
that the A— 8 Be r.m.s. distance in the \Be ground state agrees with that between the A 
and the cm. of the two en's in the aaA model calculation for A Be (2.85 fm vs. 2.88 fm, 
respectively). This three-body calculation for ^Be, as is evident from the values given 
in parentheses in Table |V11| , gives smaller Baa values by about 1.7 MeV and also smaller 
geometrical sizes than those given by the four-body calculation. The difference is due to 
the aa correlations which are absent in the three-body calculation, and which are built in 
within the Yakubovsky equations of the four-body calculation. By breaking up 8 Be into two 
a's in the four-body calculation, substantial attraction is gained due to several additional 
bound subsystems such as for the AA-He - a and the A He - A He clusters which almost saturate 
the corresponding rearrangement channels (a A A) - a and (<xA) - (<xA), respectively. The 
prominence of these channels is exhibited in Table [TX] where the probability weights of 
the various rearrangement-channel components of the AA-^e ground-state wavefunction are 
listed, using the generic expression (^). It is striking how little the (aa) - (AA) rearrangement 
channel contributes to the four-body calculation. Other calculations J^|| which pointed out 
the linear relationship discussed above did not account for the full range of 'attractive' 
rearrangement channels and consequently found smaller values, not exceeding 0.5 MeV, for 
the AB AA (^Be) gain as Vaa — ► 0. 



B. The A = 5 AA hypernuclei 

The A = 5 isodoublet hypernuclei Aj \H and AA^e were considered by us as three-cluster 
systems 3 H AA and 3 He AA, respectively. The results of Faddeev calculations for these 
species were included in Table [VII]] , in terms of the incremental binding energy AI?aa (p6|), 
together with results for A AHe and AA Be. We note that AB\a for ^jjHe is slightly larger than 
for aa-H, owing to the tighter binding of the core A hypernucleus A He as compared to \H. 
The A = 5 AA hypernuclei are found to be particle stable for all the AA attractive potentials 
used in the present calculation. As the AA interaction strength goes to zero, AB^(A = 5) 
approaches approximately zero, similarly to AB^(A = 6), provided A£?aa(^4=5) is defined 
with respect to the (2J+1) average of the ( A H, A He) levels of Table fVJ (this is a physically 
meaningful definition of AB\\ when several spin-flip nuclear core levels contribute to the 
Faddeev calculation; see also the comment made following Eq. ([/]) for the wavefunction 
of these A = 5 AA hypernuclei). We recall that both of the A = 5, 6 AA hypernuclear 
systems are treated here as three-body clusters. A nearly linear correlation appears between 
AB\^(A = 6) and A£>aa(^4 = 5), as shown in Fig. |J Judging by the slope of the straight 
lines prevailing over most of the interval shown in the figure, the A A interaction is more 
effective in binding AA He than binding either one of the A = 5 AA hypernuclei. For this range 
of mass values, the heavier the nuclear core is - the larger AB^a is (except for a vanishingly 
weak Vaa where the trend is opposite), implying that no saturation is yet reached. Within 
such three-body models A£?aa appears to get saturated already about A = 10, judging 
by the slope of the three-body calculation straight line of Fig. [3]. However, the special 
four-body cluster structure of ^Be reverses this trend of saturation, apparently delaying 
it into somewhat heavier species (for considerably heavier nuclear cores AB^a should start 
occasionally decreasing slowly to zero with A). 

We have already compared our calculations, discussing Tables [V| JVT| , |VTI| , with those of 
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Ref. [14]]. In Table [VIII] we also added a comparison with the calculations of Ref. |T5 . 
Their AI?aa values are incredibly higher everywhere than ours, and this holds also for their 
ND-type calculation of |He. 



IV. SUMMARY 

In this work we studied light AA hypernuclear systems which may be described in terms of 
few-cluster systems and treated by solving the three-body Faddeev and four-body Faddeev- 
Yakubovsky s-wave differential equations. For aa-Bg, the Faddeev- Yakubovsky solution of 
the aaAA problem given here is the first one of its kind. We estimate an error of about 0.5 
MeV due to limiting the AA-^e calculation to s waves. Our calculation yields more binding 
for AA-Be than that produced by most of the previous calculations, wherever comparison is 
meaningful. In particular, a fairly large value of AB\\ ~ 1.5 MeV survives in the limit 
Vaa - > owing to the «a structure of the 8 Be core. Our calculation confirms, if not aggra- 
vates, the incompatibility of the 'old' experimental determination of the binding energy of 
A^He [0] with that of AA^e [f2| in accordance with previous cluster calculations, irrespective of 
whether or not three-body Acta interactions are included (see e.g. Refs. [f7|,14| respectively). 
The 'new' experimental determination of the binding energy of AA-He 0] is found to be even 
more incompatible with that of AA-^e. Assuming that the determination of B^(^Be) was 
plagued by an unobserved 7 deexcitation involving either AA-^e* or A-^e*, somewhere along 



the sequence of tracks observed in the emulsion event carefully reanalyzed in Ref. JT7| , does 
not satisfactorily resolve this incompatibility. Adding A aB || as input does not alleviate it 
either, since the possibility of unobserved 7 deexcitation cannot be dismissed also for this 
species, while on the theoretical side the analysis of A aB in terms of a few-cluster system is 
more dubious than for the lighter AA species. Our conclusion differs radically from the re- 



cent claim of compatibility by Albertus et al. [40] who used a different methodology, largely 



untested in studies of very light nuclear species. 

Discarding the past history of emulsion experimentation for AA hypernuclear events 
identified as heavier than ^j(Re, because of the ambiguities mentioned here, one remains with 
the very recent report from the KEK E373 experiment [0] which claims to have identified 
uniquely AA-He, with A£>aa ~ 1 MeV. No particle-stable excited states are possible for this 
species or for its A hypernuclear core \H.e, so this event - if confirmed by adding more 
events of its kind - should be taken as the most directly relevant constraint on the AA 
interaction. Moreover, A ®H.e is ideally suited for three-body cluster calculations such as the 
Faddeev equations here solved for the aAA system. Using s-wave soft-core AA potentials 
that simulate several of the Nijmegen AA interaction models, we have shown that such 
simulations of model NSC97, versions e and / of which have been shown recently to agree 



quantitatively with data on light A hypernuclei |18 , 4l| , are capable of nearly reproducing 



the new value || of AB^( A ^B.e), short by about 0.5 MeV of the mean value.f] In fact, we 



*in order to precisely reproduce this mean value, the parameter 7 needs to be increased from its 
values listed in Table |l| for model NSC97 to 0.6598, yielding the following low-energy parameters: 
aAA = —0.77 fm, taa = 6.59 fm. 
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estimate the theoretical uncertainty of our Faddeev calculation for AA He as a fraction of 0.5 
MeV. Two approximations give rise to this theoretical uncertainty, one is the restriction to 
s-waves in the partial-wave expansion of the Faddeev equations, excluding higher i values; 
the other one is ignoring the off-diagonal AA — EN interaction which admixes E components 
into the ^He wavefunction. The following estimates may be given for the errors incurred in 



these approximations. As for the s-wave approximation, the calculations by Ikeda et al. f42 
give 99.7% probability weight for the (£, A) = (0,0) channel, with higher partial waves adding 
only 0.3 MeV binding energy, well within the error bars of the measured energies. As for 
including explicitly the AA — EN coupling, a recent work by Yamada and Nakamoto [TBI 
using model ND(G) (another model) finds an increase of 0.4 (0.1) MeV in the calculated 
_B AA ( AA He) value due to including explicitly the AA — EN coupling, resulting in a 0.3% 
(0.1%) probability E component. Carr et al. [f|3] found a similar increase of 0.2 MeV due 
to explicitly including this coupling, for a AA interaction comparable to those used in the 
present work. However, these authors pointed out that, with respect to a calculation using 
an effective single-channel AA interaction which implicitly includes the effect of AA — EN 
coupling in free space, an explicit two-channel calculation for AA He yields a reduction of 
-^AA( AA He) by about 0.3 MeV. This reduction is expected due to the tight binding of the a 
cluster which inhibits the effectiveness of the AA — EN coupling. Therefore, the effects of 
the two approximations here discussed tend to cancel each other almost completely. This 
conclusion should also hold qualitatively for AA Be due to the fairly large energy which the 
A A — EN coupling requires in order to break an a constituent. 

If model NSC97, as our few-body calculations suggest, provides for the right SU(3) 
extrapolation from fits to NN and YN data, then it is questionable whether AA H is particle 
stable. The subtlties involved in estimating whether or not this species provides for the 
onset of AA binding in nuclei were clearly demonstrated in Ref. MM. Its AB AA , at any 



rate, should not exceed 0.5 MeV for considerably stronger AA interactions | 15[ . If the AGS 
E906 events conjectured in Ref. |TU| as evidence for AA H are confirmed as such in a future 
extension of this experiment, this four-body system pnAA would play as a fundamental role 
for studying theoretically the hyperon-hyperon forces as A H (pnA) has played for studying 
theoretically the hyperon-nucleon forces ^| . Our preliminary results, studying the solutions 
of the Faddeev-Yakubovsky equations for the pnAA system, indicate that AA H is unstable 
against A emission. 
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TABLES 



TABLE I. Parameters of the potential (24) simulating AA 1 Sq Nijmegen interaction potentials. 



A (fm) 



,(*) 



(MeV) 



1 

2 
3 



1.342 
0.777 
0.350 



-21.49 
-379.1 7 
9324 



TABLE II. AA 1 So scattering lengths a and effective ranges r for several Nijmegen potential 
models. The values in parentheses correspond to using the form (^4|), with parameters from Table 
| plus the values of 7 listed here to simulate these potentials. 

Parameter NSC97b NSC97e NSC97f ND ESC00 

oaa (fm) -0.38 (-0.38) -0.50 (-0.50) -0.35 (-0.36) -2.80 (-2.81) -10.6 (-10.6) 
r AA (fm) 10.24 (15.2) 9.11 (10.6) 14.68 (16.6) 2.81 (2.95) 2.7 (2.23) 

7 0.4804 0.5463 0.4672 1 1.2044 



TABLE III. Parameters of the Isle potential ( p5| ) for the aa and Aa interactions. 



V4 J P J (MeV) (fm) Vg> (MeV) ^ (fm) 



aa [29] 


120.0 


1.53 


30.18 


2.85 


Aa |3C] 


450.4 


1.25 


404.9 


1.41 



TABLE IV. Binding energies -Ba of A hypernuclei. 







A H(0+) 


iH(l+) 


iHe(0+) 


iHe(l+) 


iHe 


A Be 


Prep 


(fm) 


1.2573 


1.2720 


1.2532 


1.2687 


1.25 




Tjcalc. 


(MeV) 


2.06 


1.04 


2.36 


1.24 


3.09 


6.64 


R exp. 
a k 


(MeV) 


2.04±0.04 a 


1.00±0.04 f ' 


2.39±0.03 a 


1.24±0.04 b 


3.12±0.02 a 


6.71±0.04 a 



a Ref. m fe Ref. m 
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TABLE V. Calculated ground-state energies of AA hypernuclei £aa (i n MeV) with respect 
to the breakup threshold of the free clusters, using aa and Aa potentials of the form (^) with 
parameters listed in Table HJ, for AA potentials (p4|) simulated by fitting to the scattering length 
a and effective range r (in fm) of several Nijmegen interaction models. 



Model 


IHe 


A Be 


OAA 


rAA 


AAHe 


AA^ e 


ESCOO 


-3.09 


-6.55 


-10.6 


2.23 


-10.7 


-19.4 


ND(G) a 






-5.37 


2.40 


-10.1 


-18.7 


ND 






-2.81 


2.95 


-9.10 


-17.7 








-0.77 


2.92 


-7.70 


-16.4 








-0.31 


3.12 


-6.98 


-15.6 


NSC97e 






-0.50 


10.6 


-6.82 


-15.4 


NSC97b 






-0.38 


15.2 


-6.60 


-15.2 


V AA = 






0.0 




-6.27 


-14.8 


ND 


-3.12 


-6.67 


-2.80 


2.81 


-9.34 


-17.15 


exp. 


-3.12±0.02 fe 


-6.62±0.04 b 






-10.9±0.6 C 
-7.25±0.19tg;:tf e 


-17.6±0.4 d 
-14.5i0.4-f 



a Ref. b Ref. Q c Ref. § d Ref. e Ref. M ^assuming A ^Be -> tt 



+ p + \Be* 



TABLE VI. Ground-state energies Eaa of ^He calculated for an aAA system; the r.m.s dis- 
tance between the A hyperons (Raa), between the center of mass of the hyperon pair to the a 
cluster (.R(aa)q) and between one of the A hyperons to the a cluster (i?Aa); the probability weight 
(P(Ao)-a) using Eq. (|j) for the rearrangement channel (Aa) — A, for various AA interaction models. 



Model 


^AA (MeV) 


Raa (fm) 


R(AA)a (fm) 


R Aa (fm) 


P(Aa)-A 


ESCOO 


-10.7 


3.09 


2.04 


2.56 


0.804 


ND(G) 


-10.1 


3.18 


2.06 


2.61 


0.830 


ND 


-9.10 


3.36 


2.11 


2.70 


0.870 


NSC97e 


-6.82 


3.93 


2.29 


3.02 


0.969 


Vaa = 


-6.27 


4.09 


2.35 


3.12 


0.999 


M ND(G) 


-9.23 


3.31 


2.14 






ND(G) 


-9.4 










PH ND 


-9.34 


3.20 




2.55 
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TABLE VII. Ground-state energies £aa of ^Be calculated for an aaAA system; the r.m.s 
distance between the a clusters (R aa ), between the A hyperons (-Raa), between one of the A 
hyperons and the cm. of the aa clusters (R( a a)A) an d between the cm. of the aa clusters and the 
cm. of the A hyperons (R( a a)(AA))i using various AA interaction models. Results for a three-body 
8 Be AA model are given in parentheses. 



Model 


E AA (MeV) 


R aa (fm) 


Raa (fm) 


R(aa)A (f m ) 


R(aa)(AA) (f m ) 


ESCOO 


-19.4 


3.3 


3.7 


2.87 


2.2 




(-17.5) 




(3.2) 


(2.49) 


(1.9) 


ND(G) 


-18.7 


3.3 


3.8 


2.91 


2.3 




(-17.0) 




(3.3) 


(2.52) 


(1.9) 


ND 


-17.7 


3.4 


3.9 


3.02 


2.3 




(-16.0) 




(3.4) 


(2.60) 


(2.0) 


NSC97e 


-15.4 


3.5 


4.2 


3.17 


2.4 




(-13.8) 




(3.7) 


(2.80) 


(2.1) 


^AA = 


-14.8 


3.6 


4.3 


3.26 


2.4 




(-13.3) 




(3.8) 


(2.85) 


(2.1) 


ND 


-17.15 


3.40 


3.02 




1.90 


© ND(G) 


(-17.7) 




(2.81) 




(1.67) 


@ ND(G) 


(-17.0) 




(2.8) 







TABLE VIII. Incremental binding energies AI?aa (hi MeV) for AA potentials fl24| ) specified by 
the scattering length a and effective range r (in fm). AB\\(A=5) is relative to the (2J+1) average 
of the ( A H, A He) levels of Table [VJ 



Model 


a AA 


fAA 


AA n 


AlHe 


A^He 


10 Bp 
AA e 


ESCOO 


-10.6 


2.23 


3.46 


3.68 


4.51 


6.1 


ND 


-2.81 


2.95 


2.11 


2.27 


2.91 


4.4 




-0.77 


2.92 


1.02 


1.13 


1.51 


3.1 




-0.31 


3.12 


0.53 


0.59 


0.79 


2.3 


NSC97e 


-0.50 


10.6 


0.50 


0.55 


0.63 


2.1 


NSC97b 


-0.38 


15.2 


0.37 


0.40 


0.41 


1.9 


^AA = 







0.11 


0.12 


0.08 


1.5 


ND 


-2.80 


2.81 






3.10 


3.74 


@ ND 


-2.80 


2.81 


2.8 


2.7 


4.3 




exp. 










4.7±0.6° 


4.3±0.4 6 



1.01±0.20±g;if c 



1.2±0.4 d 



a Ref. § 6 Ref. [0 c Ref. 



assuming AA Be — > 7r + p + 



A Be* 
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TABLE IX. Probability weights of rearrangement channels for /^Be, evaluated using a gen- 
eralization of Eq. within a four-body aoAA calculation. In parentheses, the weights for the 
bound-system clusters are given: ^A-He for aAA, ^Be for aaA and ^He for a A. 



Model P(aAA)~a P(aaA)-A P(aA)-(aA) P(aa)-{AA) 

ESCOO 0.450 0.314 0.227 0.0095 

(0.448) (0.312) (0.224) 

ND 0.379 0.359 0.253 0.0083 

(0.378) (0.357) (0.250) 

NSC97e 0.276 0.428 0.291 0.0048 

(0.276) (0.424) (0.287) 

V\a = 0.260 0.445 0.290 0.0046 

(0.247) (0.442) (0.286) 
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FIG. 1. AA-He and /^Be ground-state energies Eg calculated by solving the s-wave Faddeev 
equations for aAA, and the Yakubovsky equations for aaAA, using model ESCOO, as function of 
the number of terms N retained in the expansion (22). A cutoff radius -Rcutoff = 19 f m was used 
in both calculations. 
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FIG. 2. AA potentials for various simulations of Nijmegen models. 
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FIG. 3. Calculated binding energies -£>aa f° r AA^e m a three-body aAA model, and for ^Be 
in a four-body aaAA model and in a three-body 8 Be AA model. 
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FIG. 4. Incremental binding energies A£>aa calculated for A = 5,6 in a three-body model. 
AB^(A=5) is relative to the (2J+1) average of the (jji, A^e) levels of Table [V. 
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